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Abstract 

We study moduli stabilization in the context of M-theory on compact manifolds with 
G2 holonomy, using superpotentials from flux and membrane instantons, and recent 
results for the Kahler potential of such models. The existence of minima with nega- 
tive cosmological constant, stabilizing all moduli, is established. While most of these 
minima preserve supersymmetry, we also find examples with broken supersymmetry. 
Supersymmetric vacua with vanishing cosmological constant can also be obtained after 
a suitable tuning of parameters. 
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1 Introduction 



The stabilization of moduli is a long-standing problem in string theory. Traditionally, at- 
tempts to solve it have included non-perturbative effects, such as gaugino condensation, that 
generate a (super)-potential for the moduli fields. There is a substantial body of literature on 
this subject, particularly in the context of heterotic string theory, see for example [1]-[18]. It 
is fairly difficult to find models with proper minima in this way, and it is probably fair to say 
that successful models constructed along these lines usually require special parameter choices 
and some degree of tuning [6] . The resulting minima are quite shallow and only separated 
from runaway directions by a small barrier [19, 20]. These problems can be traced back to 
the nature of non-perturbative superpotentials, which are double-exponential in (canonically 
normalized) moduli fields. 

Recently, progress has been made by using flux of anti-symmetric tensor fields to stabilize 
moduli [21]. Unlike non-perturbative superpotentials, superpotentials from flux are merely 
single-exponential in canonically normalized moduli fields and are, therefore, more likely to 
produce minima. It has been known for some time [3, 22, 23, 24, 25] that a combination 
of flux and non-perturbative effects can stabilize moduli successfully, and under relatively 
generic conditions, but only with recent advances in the understanding of flux compactifica- 
tions [26]-[38] has this possibility been analyzed more systematically [39]-[41]. In this paper 
we will present the first systematic analysis of this problem for M-theory on manifolds of G2 
holonomy with a superpotential from flux and membrane instantons. 

M-theory on (compact) manifolds X of G2 holonomy leads to non-chiral [42] four- dimensional 
effective theories with N = 1 supersymmetry [43, 44]. Non-Abelian gauge groups and matter 
fields, which may account for the standard model particles, arise when the G2 space develops 
singularities [45]- [53]. In this paper we will not consider such possible matter field sectors, 
but focus on the gravity /mo dull part of the theory, which contains 6^(X) chiral moduli mul- 
tiplets . Their real parts parameterize the moduli space of the G2 manifold, while their 
imaginary parts corresponds to axions. We are interested in stabilizing these moduli fields 

combining the effects of flux and membrane instantons. The flux superpotential has been 
computed in Ref. [27, 44], while the result for the structure of membrane instanton contribu- 
tions to the superpotential can be found in Ref. [54]. The scalar potential of four- dimensional 

= 1 supergravity also depends on the Kahler potential for which we rely on the results 
of Ref. [55]. We will be focusing on compact G2 manifolds constructed by blowing up the 
singularities of G2 orbifolds [56, 57, 58, 45]. For such G2 manifolds, the moduli naturally 
split into two classes, {^^) = (T"',W), namely the "bulk" moduli associated with the 
underlying torus and the "blow-up" moduli W that arise from blowing up the singularities. 
We will address the stabilization of both types of moduli for realistic models. 

A practical problem is that there are no simple compact G2 manifolds with, say, 6^(X) = 1 
available. In particular, the calculation of the Kahler potential in Ref. [55] has been carried 
out for an example with 6^(X) = 43. We will, therefore, have to deal with a large number of 
moduli. This task will be approached by starting with relatively simple, but characteristic, 
toy models and then working our way up to include the full complications of more realistic 
cases. After setting up the structure of the models in the next section. Section 3 discusses 
the stabilization of the bulk moduli T"", first for a simple universal case, and then including 
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all bulk moduli. Section 4 includes the blow-up moduli W, again starting with a simple 
universal toy model and then moving on to include all moduli. Conclusions and further 
directions are presented in Section 5. 

2 General structure of four-dimensional effective the- 
ories 

In this section, we will review the structure of four-dimensional = 1 theories, which orig- 
inate from M-theory on (compact) G2 manifolds. These four-dimensional models will form 
the basis for the subsequent analysis in this paper. For mathematical facts on manifolds with 
G2 holonomy we refer to Ref. [58]. Many of the physics aspects of the following presentation 
can be found in Refs. [43, 44]. 

Before we specialize to the models considered in this paper, let us collect a number of 
useful general facts about M-theory on G2 manifolds. A seven-dimensional manifold X with 
holonomy G2 has a vanishing first Betti number, b^{X) = 0, and its cohomology is, therefore, 
characterized by b'^{X) and b^{X). When compactifying M-theory on such a manifold, the 
resulting four- dimensional low-energy theory has = 1 supersymmetry and contains b'^{X) 
Abelian vector multiplets and b^{X) uncharged chiral multiplets [42, 43]. We denote these 
latter chiral fields, which are the primary focus of this paper, by I, J, ... = 1, . . . , b^{X). 

To compute the scalar potential we will need an explicit expression for the Kahler poten- 
tial of these chiral superfields. Ideally, one would like a generic and explicit formulae that 
may depend on topological data of the G2 manifold in question, but still applies to all G2 
manifolds. Unfortunately, no such formula is known. We will go around this problem by 
focusing on a particular class of G2 manifold for which an (approximate) form of the Kahler 
potential is known. However, let us first collect some general facts. It is known [44] that the 
Kahler potential can be expressed in terms of the volume V of the internal manifold as 

ir = -31n(^). (2.1) 

This relation implies that K only depends on geometrical data and is, hence, independent 
of the axionic fields, which we take to be the imaginary parts of the chiral fields We, 
therefore, have 

K = K{^^ + ¥) . (2.2) 

Further, it can be shown [44] that a rescaling of the G2 space by a factor A, corresponding 
to a change V — » A^V of the volume, leads to a scaling + A^($'^ + $'^). Hence, the 

real parts of the chiral fields $^ measure three-dimensional volumes. More precisely, they 
are proportional to the volumes of (a basis of) three cycles of the G2 manifold. 

It can be shown by straightforward dimensional reduction [55, 43] that the gauge kinetic 
functions / of the gauge multiplets are of the form 

/ = c,$' , (2.3) 

where C/ are constant coefficients that depend on the G2 manifold and the gauge multiplet 
under consideration. 
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Perturbatively, and in the absence of flux, the superpotential W vanishes. It has been 
shown [27, 44] that internal flux of the M-theory four-form field strength leads to a super- 
potential 

Wfiux = mi^' , (2.4) 

where the flux parameters mj are quantized. Non-perturbatively, an additional contribution 
to the superpotential is generated by membrane instantons [54] wrapping three-cycles within 
the G2 manifold. A typical instanton superpotential is given by a sum of terms of the form 

VTinst ~ exp(-pj$^) , (2.5) 

where the constants pi determine the three-cycle in question. In this paper we will study 
the stabilization of G2 moduli due to a combination of flux and membrane instanton con- 
tributions to the superpotential. Note that the flux terms (2.4) in W break the discrete 
shift symmetries of the axions that, even after the inclusion of instanton effects, one would 
normally expect. Hence field values that differ by integer multiples of 27r2 are no longer 
equivalent, as would have been the case in the absence of flux. In Ref. [44], these integer 
shifts of the axions have been related to the possible values h of the external part of the flux. 
More precisely, the external flux is related to the axionic components of the superfields by 

h ~ V-^mj5>(<l>^) , (2.6) 

where we have dropped an overall constant. 

We will focus on a particular class of compact G2 manifolds, as constructed by Joyce [56, 
57, 58]. These G2 manifolds are built from seven dimensional G2 orbifolds with co-dimensional 
four singularities by blowing up their singularities. Accordingly, the moduli $^ split into two 
classes, namely the "bulk moduli" T", where a,b, . . . = 1, . . . , A, associated with the under- 
lying orbifold and the "blow-up moduli" W, where = 1,...,/. Geometrically, the 
real parts of the bulk moduli measure the radii (or angles) of the underlying torus, while the 
blow-up moduli measure the size and orientation of the blow-ups. It will be convenient to 
split the moduli fields into their real and imaginary parts by writing 

= r + ir" , W =u' + iu' . (2.7) 

Recall that, more precisely, the real parts and are proportional to volumes of internal 
three cycles, while the fields and i/* are axions. We normalize the real parts so they 
measure these three-cycle volumes in units of T2~^, the inverse membrane tension. 

Here we consider examples where the only moduli of the underlying orbifold are its seven 
radii. Hence, we have seven bulk moduli, each with a real part proportional to the volume 
of one of the seven bulk three-cycles. The calculation of Ref. [55] shows that the Kahler 
potential for these seven bulk moduli is given by 

7 

i^buik = -5^1n(T'^ + f'^) , (2.8) 

a=l 

as expected. The idea is to obtain an approximate Kahler potential for all moduli by ex- 
panding in the typical size p of the blow-ups. Blowing up co-dimension four singularities 
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introduces additional two-cycles, and the blow-up moduli should, therefore, scale with 
p^. The construction of Ref. [56, 58] also shows that the leading correction to the orbifold 
volume due to the blow-ups arises at order p^. The leading correction to Eq. (2.8) must, 
therefore, be proportional to {W + f/*)^. The general Kahler potential can then be written 
as [59] 

7 / 

= -^ln(T" + T") +^/i(f/* + f7*)2 , (2.9) 

a=l i=l 

where /, are functions of the real parts of the bulk moduli. Under a rescaling of all 
moduli with the first term in Eq. (2.9) already reproduces the correct scaling of K, or 
equivalently, via Eq. (2.1), of the volume V. Consequently, the second term in Eq. (2.9) 
should be invariant under the rescaling. This implies the functions must be homogeneous 
of degree —2 in t". For practical purposes, we will use the simple form 

/, = 8[nLi(T'^ + TTi"' , (2.10) 
where the vectors Pj = (pia) satisfy 



a=l 



P^a = 2, Pia e {0,1} (2.11) 



for all i. This form has indeed been found in the explicit calculation of Ref. [55] and can 
be expected to hold for a wider class of examples. We will use the Kahler potential (2.9), 
with the functions /j given by Eq. (2.10), as the basis of our analysis. We should stress 
that this Kahler potential is approximate and only holds for sufficiently small blow-ups. 
More precisely, one expects corrections of order u^/f^ to Eq. (2.9) and one should, therefore, 
require the ratio u/t to be smaller than one. We will discuss the constraints on the validity 
of Eq. (2.9) in more detail below. 

In this paper we keep the number / of blow-up moduli arbitrary. For the examples 
presented in Ref. [58] it ranges from a few to well over 100, but we will see later that 
consistent minima can sometimes only be obtained for sufficiently low I. The structure 
of the powers pia will not be of particular importance for our analysis and, apart from 
the constraints (2.11), we will keep them arbitrary. Despite this level of generality, it will 
be useful for us to present a concrete example for illustration. This example is based on 
the orbifold T^/Zg, and contains twelve co-dimension four singularities locally modeled on 

X C^/Z^. For the details of the construction we refer to Refs. [55, 56, 58]. The singularities 
are removed by blowing-up the origin in C^/Z^. Ignoring some subtleties associated with 
the smoothing required to join them to the rest of the manifold (see [55] for details), the 
blow-ups support an Eguchi-Hanson metric 



^•seh — :: ;x + 



1 — ^ 



r 
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1 7 ] + ao + a. 



j,4 



where a* are Maurer-Cartan one forms, and the parameter p is the radius of the blow-up. 
Each blow-up is associated with three moduli (giving I = 36) whose real parts have the form 
3?(?7*) ~ p^(n")^, where n is a unit three-vector describing the blow-up's orientation. The 
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P(r,a) 


a = 1 


a = 2 


a = 3 


r = 1 


(1,0,0,0,0,1,0) 


(0,1,0,0,1,0,0) 


(0,0,1,1,0,0,0) 


r = 2 


(1,0,0,0,0,0,1) 


(0,0,1,0,1,0,0) 


(0,1,0,1,0,0,0) 


r = 3 


(1,0,0,1,0,0,0) 


(0,0,1,0,0,1,0) 


(0,1,0,0,0,0,1) 



Table 1: Values of vectors P{r,a) that define the moduli Kahler potential for the G2 manifold 
based on T'^/Zg. 



real parts of the seven bulk moduli, on the other hand, are products of three of the radii of 
the underlying torus, chosen so that they are the volumes of the three cycles preserved by 
the orbifolding. The 36 blow-up moduli can be labeled by a triple of indices (i) = (r, n, a) , 
where r = 1, 2, 3, n = 1, 2, 3, 4 and a = 1, 2, 3. The index r refers to the three Z2 symmetries 
of the orbifold and it indicates that a given blow-up modulus is associated with a fix point 
under the t^^ Z2 symmetry. Further, a is an orientation index and n numbers the moduli 
of a given fixed point type and given orientation. The powers pj only depend on the type r 
and the orientation a and the nine relevant vectors P(T,a) are given in Table 1. 

We will consider a superpotential W that arises from a combination of flux and membrane 
instantons. From the above discussion this should have the form 

7 I 
W = J2 ^i^a, K: T'^) + ^(Z'^' ' (2.13) 

a=l 1=1 

where the function w is defined by 

w{m,k,X) = mX + ke~^ . (2.14) 

The parameters rria and fii are the flux parameters for the bulk and blow-up moduli, respec- 
tively. They are quantized and we will take them to be integers for simplicity. This can 
always be arranged by absorbing the overall factor into the normalization of the potential 
and rescaling ki, U. Strictly, the pre-factors ka and U for the instanton contributions are 
(non-exponential) functions of the moduli. However, these functions have not been explic- 
itly computed for the cases at hand and we will take them to be constant and real, if only 
to make the problem manageable. 

To summarize, the models we are going to consider contain seven bulk moduli and an 
arbitrary number / of blow-up moduli, with a Kahler potential (2.9)-(2.11) and a superpo- 
tential (2.13), (2.14). Let us now discuss the constraint we should impose on moduli space. 
For the supergravity approximation, which underlies the above expressions for K and W , 
to be valid we need all internal length scales to be larger than the 1 1-dimensional Planck 
length. Given the normalization of our moduli this means that we should require 

> 1 , u'>l . (2.15) 

In addition, to suppress higher 0{uf) terms in the Kahler potential, we should also ensure 
that 

4$^/.wN<l- (2.16) 

i=l 
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Geometrically, this constraint implies that the fraction of the orbifold volume "taken away" 
by the blow-ups is small or, in other words, that we are not too far away from the orbifold 
limit. While we clearly need to impose a constraint of this type it is, in practice, difficult to 
decide how small the left-hand side of Eq. (2.16) really has to be for the unknown corrections 
to the Kahler potential to be negligible. We will later adopt a more rigorous approach 
whereby we add hypothetical 0{uf) terms to the Kahler potential and check whether they 
leave the lowest order results essentially unchanged. We will then only accept results that 
pass this consistency test. Note that Eq. (2.15) provides a lower bound on individual while 
the left-hand side of the constraint (2.16) scales with their total number I. We, therefore, 
anticipate that consistent results are more difficult to obtain as the number of blow-up moduli 
increases. 

In this paper we are interested in the scalar potential for the fields ($'^) = (T", [/*), which 
is given by the standard = 1 supergravity expression 



(2.17) 



where Kjj is the Kahler metric 

"■''WW <^-'«> 

and K^"^ its inverse. The F-terms Fj are defined by 

Fi = Wi + KiW , (2.19) 

where indices on W or K generally indicate derivatives, that is, 

dK dW 

Let us summarize the main general properties of this = 1 supergravity potential, relevant 
to our applications. A minimum of the potential only respects supersymmetry if F/ = for 
all /. Conversely, a solution to the F-equations 

Fi = (2.21) 

constitutes a supersymmetric extremum of the potential. The second derivatives of at a 
supersymmetric extremum are given by 

Vij\p^=o = -e'^WFjj (2.22) 



K'''^FkiFij-2Kjj\W\'' . (2.23) 



In general each individual case has to be checked to see if this matrix is positive definite and 
such an extremum is indeed a minimum. Suppose that we have a family of potentials that 
depend on a set of parameters and that a minimum has been found for a particular point in 
parameter space. The fact that the Jacobi matrix of V is positive definite, and hence non- 
singular, at this minimum, together with the implicit function theorem, then tells us that 
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minima exit in an open neighbourhood of this particular point in parameter space. In some 
of our subsequent examples, we will explicitly find minima for special non-generic values of 
measure zero in parameter space. The above statement then establishes the existence of 
minima for more generic parameter choices in an open neighbourhood of parameter space. 
The potential value at a supersymmetric extremum is 



and, hence, the cosmological constant is either negative if {W ^ 0) or zero if {W = 0), 
leading to AdS4 or four- dimensional Minkowski space, respectively. Conversely, a positive 
cosmological constant implies broken supersymmetry. Supersymmetric extrema with van- 
ishing cosmological constant are characterized by 



and are, hence, independent of the Kahler potential. They are only minima if Wjj is non- 
singular, as can be seen from Eqs. (2.22), (2.23) by setting W = 0. Again, suppose we 
have a family of superpotentials depending on a set of parameters. Finding a solution to 
the equations (2.25) will usually imply a fine-tuning in parameters space, as is normally 
required to set the cosmological constant to zero. Now consider the first derivatives Fjj and 
Fjj at a solution of (2.25). They are given by Fjj = Wjj and Fjj = and, hence, the 
block matrix consisting of Fjj, Fjj and their conjugates is non-singular. From the implicit 
function theorem this establishes the existence of supersymmetric minima (generically with 
negative cosmological constants) in an open neighbourhood in parameter space, around the 
special point where the cosmological constant vanishes. 

3 Stabilizing bulk moduli 

The structure of our models is fairly complicated and finding minima of the potential in 
the general case would be a quixotic task. The strategy we will adopt starts with simplified 
models, gradually working our way up to include the full complexity of the system. As we 
will see, the results for the simplified models can be used to construct minima in the general 
case. In this section, we focus on the bulk moduli T° and neglect the blow-up moduli U\ 
They will be included in the subsequent section. 

3.1 The universal case 

As a further simplification, we first consider a "universal" case where we set all seven bulk 
moduli equal. Such a model then contains a single superfield 




(2.24) 



w 



0, 



(2.25) 



T = t + iT . 



(3.26) 



with Kahler potential 



K 



7in(r + r) 



(3.27) 



and superpotential 





7 



We can think of t as the breathing mode of the G2 manifold. Of course, we are not claiming 
the above model corresponds to a particular G2 manifold in any strict sense. It is merely a 
toy model that incorporates the main features of the bulk sector of our realistic models. In 
addition, we will see that the results from this model can be transferred to realistic cases. 

With this understood, we are now going to analyze the properties of the above model. 
The F-term can be written as 



7 

w'{m, k, T) — —w{m, k, T) 



(3.29) 



where w' is the derivative of w with respect to its last argument. More explicitly, split up 
into real and imaginary part, this reads 



3?(Fr) 
3(Ft) 



m 
' 2" 



K 2 



7 



cos r e 



where we have defined 



m 








K (^24 


■D 


2" 










k 




K = 





sm r e 



7r' 

T 



(3.30) 
(3.31) 



m 



(3.32) 



The scalar potential (2.17) for the Kahler potential (3.27) and a general W is given by 



7. 



7 



V = — ■> iWrr - -^{WtW) + -\W 
while inserting the explicit form (3.28) of W leads to 



C2 



7r2 



5 cosfr) 



-K 



■ 7 7' 



-2t 



. , .J 14rsin(r) 
[cos(r) — rsm(rj) 



t 



(3.33) 



(3.34) 



Here the overall constant is given by 



m 
224 



(3.35) 



It can be shown from the F-terms (3.30), (3.31) that every supersymmetric extremum is 
indeed a minimum of the potential. Conversely, we have found numerically that every 
minimum of the potential (3.34) preserves supersymmetry. It is, therefore, sufficient to look 
at the F-equations. 

It is easy to show from Eqs. (3.30) and (3.31) that the solutions to the F-equations can 
be characterized as follows. The values of t at the minima are the solutions to the equation 



cos/(t) 



1 5t 



(3.36) 
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where 

/(t) = ^V/t2(2t + 7)2e-2*-25t2 , (3.37) 

for which the signs of k and sin f{t) are the same. For each such t, there are two associated 
solutions for the imaginary part r, namely 

T = ±f{t) . (3.38) 

From Eq. (3.36) all (positive) solutions for t must be in the range 

te[o,w], (3.39) 



where tmax is defined by 



Very roughly we have 



5tmax (340) 



2f -\-7 



W ~ In |k| . (3.41) 

In order to have solutions that satisfy the constraint (2.15) on t, we need the parameter \k\ 
to be sufficiently large, typically of order one or larger. We note that while t varies between 
zero and tmax the function / varies between \k\ and zero. This means, from Eq. (3.36), that 
the typical number, A'soi, of minima is given by 

iVsoi — , (3.42) 

TT 

where we have remembered the two possibilities (3.38) for the imaginary part r. This 
relation should, of course, be understood as a rough indication of the number of minima for 
sufficiently large \k\. It shows that, for \k\ ^ 1, their number is substantial. For k < we 
have a particular minimum with vanishing imaginary part at 

t = tmax , r = . (3.43) 

All other minima are complex and t < tmax- These statements can be easily checked in 
Figure 1, where we show a contour plot of the scalar potential, as given by Eq. (3.34), in 
the t, T plane. We have chosen m = —1, k = 20, and we can see that, following Eq. (3.42), 
the number of minima is very close to l^l/vr, with n = k/m = —20; moreover, and according 
to Eq. (3.41), the value of tmax is given by ln|K| and, also, we can see in the figure that 
Eq. (3.43) is fulfilled. Finally we have checked that all these minima preserve supersymmetry 
with negative vacuum energy. 

More precisely, from Eq. (2.24) the potential value at any one of these minima is given 

by 

y. ,_21C! (t+l)- + r^ 

and is, hence, always strictly negative. Also, since from Eq. (3.38) r increases as t decreases, 
the absolute minimum will be the one with the smallest value of t. For k > 0, or sufficiently 
negative n ^ —5, this is always a complex minimum, r 7^ 0. Also note that the two minima 
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0.5 - 



qI 1 1 1 1 1 1 

-6 -4 -2 2 4 6 

x/k 

Figure 1: Contour plot of the potential, given by Eq. (3.34), in the (t,r) plane for m = — 1, 
k = 20. 



for a given solution of t are degenerate, since they only differ by a sign of the imaginary part 
r. 

To summarize, for our simple model, the structure of minima is controlled by the parame- 
ter K, defined in Eq. (2.9). To have minima compatible with the supergravity approximation, 
that is, minima satisfying Eq. (2.15), we need \k\ ^ 1. The number of minima is roughly 
given by l^l/vr and the real part t is always in the range t G [0,tmax], where, roughly, 
^max ~ In For K < a special real minimum with t = tmax and r = exists; all other 
minima are complex. All minima preserve supersymmetry and have a negative cosmological 
constant. The absolute minimum is the one with the lowest value of t, and generally has a 
non- vanishing imaginary part, r 7^ 0. 

3.2 Non- universal bulk moduli 

We now move on to a more realistic model with seven individual moduli 

T'^ = ta + tTa, (3.45) 

where a, fe, . . . = 1, . . . , 7. The real parts ta can be interpreted as the seven radii of the torus 
that underlies the construction of the G2 manifold. The Kahler potential and superpotential 
are now given by 

7 

K = -^ln(T'^ + f'^) (3.46) 

a=l 
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W = Y.'^ia), (3.47) 

a=l 

where the function w is defined in Eq. (2.14) and we have written 

W(a) = w{ma,ka,T°-) , (3.48) 
for ease of notation. The F-terms take the form 

1 

while the scalar potential can be written as 



Fa = w[,) - —W , (3.49) 



V = ^^1^ |e - tMn^UW)] + IW^P} • (3.50) 

Here w' denotes the derivative of w with respect to its last argument. Analyzing the general 
vacuum structure of this model would be rather cumbersome given the large parameter space. 
We will not attempt to do this but instead present a number of examples. 

Consider a "universal" choice of parameters, independent of the index a, that is, 

rUa = m , ka = k (3.51) 

for some real m, k and all a = 1, . . . , 7. Then the F-terms (3.49) can be written as 

1 7 

Fa = w'{m, k, T^-T^Yl ^' • (3.52) 



b=l 



Now take any of the supersymmetric minima T = t + ir found in the previous subsection 
for the universal case. Comparison of Eqs. (3.29) and (3.52) then shows that setting 

T" = T , (3.53) 

for all a = 1, . . . , 7 solves all seven equations Fa = 0. Hence, every supersymmetric minimum 
of the universal single T model leads to a supersymmetric extremum for the case with seven 
moduli Ta, provided the parameters have been chosen universally, as in Eq. (3.51). It can be 
shown from Eqs. (2.22) and (2.23) that these extrema are indeed minima if t is sufficiently 
large, that is, roughly t ^ 1. As in the universal case, the cosmological constant is always 
negative. As discussed at the end of Section 2, these minima still exist for a non-universal 
choice of parameters sufficiently close to the universal one. 

Allowing seven T-moduli also opens up qualitatively new possibilities compared to the 
universal case. Consider the conditions (2.25) for supersymmetric minima with vanishing 
cosmological constant. Applied to the present case they lead to moduli values 

T'' = \n — + TxiHa (3.54) 
ma 
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where Ua is an integer which is even for positive ka/ma and odd otherwise. In addition, 
vanishing of the real and imaginary part of the superpotential imphes that 



a=l 



In 



rrin 



a=l 



maUn 







. 



(3.55) 
(3.56) 



The first of these conditions is the usual fine-tuning required to make the cosmological con- 
stant vanish. It constrains parameters to a particular hyper-surface in parameter space. 
With favourable choices for the signs of ka/ma, the second condition will have an infinite 
number of solutions. They correspond to an infinite number of degenerate supersymmetric 
minima with vanishing cosmological constant that differ by integer shifts in the axion di- 
rections. We present an example of this behaviour in Figures 2,3, where we have split the 
moduli into two types, a = 1, 2, 3 and b = 4,5, 6, 7, in order to make the graphics manage- 
able. In Figure 2 we show a contour plot of the potential, as given by Eq. (3.50), in the ta, tf, 
plane, for values of the parameters rria = 4, ka = 200, rrib = —3, kb = —150, with the imag- 
inary parts of the fields, Ta, Th fixed at zero. In this way, both Eqs (3.55,3.56) are fulfilled. 
Moreover, a minimum is found at the predicted values for the moduli, given by Eq. (3.54) 
with Ua = Ub = 0, which, in this case, becomes ta = tb = 3.9. We can also see that for larger 
values of ta a maximum appears, while for larger tb we have a saddle point. In Figure 3, 
where we show the potential as a function of Ta, Tb for ta, tb fixed at their minimum values, 
we can check how Ta = Tb = does indeed correspond to a minimum which, however, is 
degenerate with those at Ta = Tb = 2mT for n integer. There exists, therefore, in these cases, 
an array of degenerate minima with zero cosmological constant. An obvious consequence of 
this pattern should be the formation of domain walls between the different minima. 

Comparison of Eqs (2.6) and (3.56) show that all these minima correspond to vanishing 
external fiux.^ From our general discussion, we also know that supersymmetric minima, 
although with negative cosmological constant, exist in a neighbourhood of the hypersur- 
face (3.55) in parameter space. 



4 Including blow-up moduli 

We should now add blow-up moduli to the models of the previous section. We begin with 
the universal case. 

4.1 The universal case 

An obvious extension of the simple model with a single bulk modulus T presented in Sec- 
tion 3.1 adds a second field U as the typical representative of a blow-up modulus. Split into 
real and imaginary parts we write 

T = t + iT (4.57) 
^We remark that minima with non- vanishing external flux can be obtained for complex parameters ka- 
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3 3.5 4 4.5 5 5.5 6 6.5 7 

Figure 2: Contour plot of the potential, given by Eq. (3.50), as a function of ta {a = 1, 2, 3) 
and tb {h = 4, 5, 6, 7) for nia = 4, ka = 200, rrib = —3, kb = —150. The imaginary parts of all 
fields have been set to zero, where we have a minimum. 



0.12^ 




-4 -4 



Figure 3: Plot of the potential, given by Eq. (3.50), as a function of r^, u for the same values 
of parameters as in the previous figure. The real parts of the fields have been set to their 
minimum values at ln|/ci/mj|, i = a,b. 



U = u + iu. (4.58) 
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The model is then defined by 



K = -71n(T + f)+6|^±££ (4.59) 
W = w{m,k,T) + w{^,l,U) , (4.60) 

where is a real positive constant. To make contact with the full model later on we will 
need to set b = 8/, where I is the number of blow-up moduli. However, for the purpose of 
this sub-section, we will treat 6 as a phenomenological parameter. It is useful to introduce 
the ratio 

e = - , (4.61) 

of blow-up and bulk modulus. Recall that the above Kahler potential should be viewed as 
an expansion in e where terms of order and higher have been neglected. In accordance 
with the general constraint (2.16) we should, therefore, work in the region of moduli space 
where 

fce^ < 1 . (4.62) 

Note that corrections of 0{e'^) to K will lead to 0{e^) terms in the scalar potential. Hence, 
on the basis of the Kahler potential (4.59), we can reliably compute the scalar potential only 
up to terms of C(e). Subsequent formulae will be quoted to this order. 

It is crucial to check that relevant features of the scalar potential, such as minima arising 
at 0{e), are stable under inclusion of terms of 0{e^) and higher, so that (4.62) is indeed 
satisfied to the necessary degree. In practice, we will add a hypothetical correction 

= (4.63) 

(T + TY ^ ' 

to the Kahler potential (4.59), where d is a real number, and compute the scalar potential 
up to order by including this correction. We will then compare the potentials at 0{e) and 
(9(e^), and only accept minima if they consistently arise at both orders and for a reasonable 
range of d. 

For the Kahler potential (4.59), and a general superpotential, the F-terms, to order e, 
are given by 

(4.64) 
(4.65) 

To this order, the scalar potential reads 

(4.66) 



Ff — 


Wt- 


2t 


Fu = 


Wu + 


t 
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The explicit expressions for the F-terms and the scalar potential are fairly complicated 
and are given in Appendix A. Inspection of the general scalar potential (2.17) shows that 



dV 



r = 0,1/ = 0) = 



dV 
dv 



r = 0,z/ = 0) = 



(4.67) 



so the potential is extremized in the imaginary directions at r = i/ = 0. Here we will focus 
on this real case, that is, we will set r = i/ = 0. Equations (A.5)-(A.8) show that then the 
imaginary parts of the F-terms are automatically zero, that is, Q'(-Ft) = '^{Fu) = 0, while 
the real parts simplify to 



3?(Ft) 



m 

T 



5 + K (^2 + e"* + ^(Mm + Le"") 



m 



be 



M - Le-'" + —{t + Ke-'] 



(4.68) 
(4.69) 



Here we have normalized our parameters with respect to the flux m by defining 



k 

K = — , 

m 



m 



m 



(4.70) 



For vanishing imaginary parts the potential (A. 10) then simplifies to 



V 



t5 



1 + K 

+7L2 



3u 
T 



1 u 1 

26^ t2 +^ 



-2m 



±LK 



7M2 u 

+ +AM- + KM 

2b t 



7 7 
Au 7 7u 14 
1v? 



~2t 



-t—u 



(4.71) 



3m 


7u 


e"* + LM 






_T 







t2 



7u 



After minimizing this potential in t and u we still, of course, have to check whether the 
extrema in the imaginary directions are indeed minima. 

Let us start by looking for supersymmetric minima with vanishing cosmological constant, 
that is, solutions to the equations Wt = 0, Wu = and W = The first two of these 
equations immediately imply that 



T 
U 



In \n\ + Tiip 



In 



L 



M 



TTiq 



(4.72) 
(4.73) 



where p (q) is an integer which is even for k > (L > 0) and odd otherwise. Vanishing of 
the superpotential leads to the two conditions 



m (1 + In \k\) + /i ( 1 + In 



mp + fiq 



0, 
. 



(4.74) 
(4.75) 
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For fixed (integer) flux m and /i the first of these equations can always be solved for ap- 
propriate choices of the parameters k and /. The second equation has an infinite number 
of other solutions when the signs of k and /x// are chosen appropriately. In particular, for 
K > and n/l > 0, we can take p = q = and obtain a real supersymmetric minimum 
with vanishing cosmological constant. Hence, by tuning in parameter space, we can have 
supersymmetric minima with vanishing cosmological constant, infinitely degenerate by inte- 
ger shifts in the axion directions. From Eqs. (2.6) and (4.75) all those minima correspond 
to vanishing external flux. In the neighbourhood in parameter space of each such minimum 
there will be supersymmetric minima with negative cosmological constant. We note that 
their existence is independent of the Kahler potential. Hence, we do not need to require 
that <^ 1, although this could be easily arranged by choosing parameters. However, 
unless be^ <^ 1, we are unable to explicitly write down the scalar potential close to these 
minima. We confirm the previous statements by showing an explicit example of a real su- 
persymmetric minimum with vanishing cosmological constant in Figures 4,5. The choice of 
parameters is such that Eq. (4.74) is satisfied, and the values of the fields correspond to 
t = ln(K) = 5.99, u = ln(//yu) = 1.1, in agreement with Eqs (4.72,4.73). In Figure 4, the con- 
ditions Ft = Fu = are also plotted to show that the minimum is indeed supersymmetric. 
Here we are considering 6 = 8, i.e. / = 1, and we have checked that corrections of order 
and higher do not affect the existence and position of the minimum. 



Let us now drop the condition of vanishing cosmological constant and focus on real 
supersymmetric minima, that is, r = = 0. The relevant F-terms have been given in 
Eqs. (4.68) and (4.69). It can be shown that solutions to the F-equations are always minima 
for the model at hand. A very rough approximation to the equation 9ft(F[/) = leads to the 
solution 

L 



u 



In 



M 



(4.76) 



which coincides with the result (4.73) one obtains for vanishing cosmological constant. Since 
we have set the phases to zero, it exists only if L and M have the same sign. Inserting this 
into the equation 1R.{Ft) = we find t is determined by 



5 + fi: 2 + 



t 



+ 



7M 



t 



In 



L 



M 



. 



(4.77) 



tmax, where 



For sufficiently small M this can be solved for negative k leading roughly to t 
tmax was defined in Eq. (3.40). Hence, we expect supersymmetric minima for k < and L, 
M having the same sign. Further, since tmax ~ In we typically need to have |k| ^ 1 in 
order to satisfy the constraint 6e ^ 1. This kind of supersymmetric minima, with negative 
cosmological constant, is shown in Figure 6 for values of the parameters given by m = — 1, 
k = 1000, /i = 3, / = 6 (i.e. L, M < 0). The minimum in the u direction is quite close to 
the value given by Eq. (4.76), whereas Eq. (4.77) constitutes a rough approximation for the 
corresponding value of t (given hj t = 8.25). Similar minima can be obtained, as already 
explained, for k < and the opposite choice of signs for L, M. And, again, we have chosen 
/ = 1 in order to keep our perturbative expansion well under control. 

For K > 0, using the relations (4.72) with p = and (4.74), which describe the super- 
symmetric minima with vanishing cosmological constant, Eq. (4.77) is identically satisfied. 
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Figure 4: Contour plot of the potential, given by Eq. (4.71), in the t, u plane, for m = 3, 
k = 1200, /i = —10, I = —30. We have also plotted the conditions Ft = Fu = in order to 
show the supersymmetric character of the minimum. The imaginary parts of all fields have 
been set to zero. 



V x10-= 



2.5^ 
2., 
1.5^ 




5.6 0.5 



Figure 5: Plot of the potential, given by Eq. (4.71), as a function of t, u for the same values 
of parameters as in the previous figure. The minimum corresponds to Eqs (4.72,4.73,4.74) 
being fulfilled. 



as it should be. Perturbing away from this special point in parameter space one can still 
find solutions to Eq.(3.36) that are characterized by k > and M < 0. They are precisely 
the supersymmetric minima close to the ones with vanishing cosmological constant whose 
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Figure 6: Contour plot of the potential, given by Eq. (4.71), in the (t, u) plane, for m = — 1, 
k = 1000, /i = 3, Z = 6. We also plot the conditions Ft = Fu = to show the supersymmetric 
character of the minimum. The imaginary parts of all fields have been set to zero. 



existence we have inferred above from general argument. 

In summary, we expect three classes of supersymmetric minima, depending on the signs 
of the various parameters but all with L/M > 0. For k < we can have either sign of M 
while for K > supersymmetric solutions only exist for M < 0. 

We now turn to the search of minima with broken supersymmetry. Unfortunately, here 
we can not count on the F-equations being fulfilled in order to look for minima, however, as 
it turns out, a good starting point is to look for solutions to the condition = in order 
to find minima of the potential (4.71). As it has already been pointed out, in order to keep 
the supergravity approximation, Eq. (2.15), valid, as well as to suppress higher order terms 
in the Kahler potential, Eq. (2.16), we need t to be large at the minimum (at least of order 
10). This translates into a large value for k. At the same time we need to keep the value of 
u of order 1, and the number of blow-up moduli below about 10, in order to comply with the 
condition be^ ^ 1. This also guarantees the smallness of the mixing terms between t and 
u in the scalar potential, and explains why the minimization along the t direction is almost 
unchanged with respect to the purely supersymmetric case. 

Therefore minima with broken supersymmetry appear for large k (and Ft ~ 0) and l,fi 
of the same order of magnitude. As for their signs, we have found only minima for / < 0, 
/i > whereas the opposite choice gives rise to a runaway potential for increasing values of 
u. A typical example of what has just been described is shown in Figures 7,8, where we plot 
both the contour and the shape of the potential given by Eq. (4.71) for m = —1, k = 1000, 
yU = l,/ = — l.Sasa function of t and u (their imaginary parts having been set to zero). 
In Figure 7 we have also plotted the constraint = in order to show how close it is 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

u 

Figure 7: Contour plot of the potential, given by Eq. (4.71), in the (t, u) plane, for m = — 1, 
k = 1000, /i = 1, / = — 1.5. We have also added the condition Ft = 0. The imaginary parts 
of all fields have been set to zero, where we have a minimum. 



V xio-= 




Figure 8: Plot of the potential, given by Eq. (4.71), as a function of t, u for the same values 
of parameters as in the previous figure. 



to the minimum (whereas Fu = lies well outside the plot). All the minima with broken 
supersymmetry that we have found have a negative cosmological constant. 

Finally we would like to add a further comment on the stability of our results. As already 
mentioned several times, we are imposing the constraint be"^ ^ 1 in order to guarantee that 
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(as yet unknown) higher order corrections to the Kahler potential will not spoil the results 
presented here. This, in turn, means that we have a tight restriction on the number of blow- 
up moduli (given by /, with b = 81) allowed in our models. We have found that, in order 
to achieve Wmin ~ 1, / = 1 becomes almost the only choice. We can still find minima with 
< Urain < 1 for / < 10, however larger values of / would result in the minima shifting to 
negative values of u. There are, though, plenty of examples in the literature of compact G2 
manifolds with a small number of blow-ups, and while the precise Kahler potential for those 
cases is not known, we expect it to resemble the one used here. In that respect our results 
should be taken from a purely phenomenological point of view. Once a complete formula for 
the Kahler potential is known, it should be very easy to incorporate it to our analysis, and 
it would be surprising if the results differed substantially from those presented here. 



4.2 The general model 

We are now ready to analyze the general model defined by the Kahler potential (2.9) and 
the superpotential (2.13), which we write as 



w 



(4.78) 



a=l 



i=l 



where W(a) = w{ma, ka, T"-) and = w(/ij, ?7*), and the function w has been defined in 
Eq. (2.14). The F-terms to order e then take the form 



= w'(. - —W 

2ta 

with /* given in Eq. (2.10). For the scalar potential we find 
1 



(4.79) 
(4.80) 
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a=l i= 
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i=l 
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(4.81) 



a=l i=l 



We will not attempt a general classification of all minima of this potential but rather 
present various classes of examples. We start with the supersymmetric minima with van- 
ishing cosmological constant, characterized by the equations w^^-^ = 0, w'^--^ = and 1^ = 0. 
The field values are easily solved for and are given by 



= In 
U' = In 



(4.82) 
(4.83) 
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where Ua are even (odd) integers for ka/ma positive (negative) and, similarly, Ui are even 
(odd) integers for li/fii positive (negative). Vanishing of the real and imaginary parts of the 
superpotential implies 



a=l ^ " ^ 1=1 ^ 

7 / 



. 



(4.84) 
(4.85) 



a=l 



i=l 



As in the analogous cases before, the first of these conditions is the usual fine-tuning in 
parameter space required to set the cosmological constant to zero, while the second has 
infinitely many solutions for favourable signs of ka/ma and h/^i (or appropriate choices of 
the flux) . We have hence found an infinite number of supersymmetric minima with vanishing 
cosmological constant that differ by integer shifts in the axion directions. As before, all those 
minima correspond to vanishing external flux, as Eq. (2.6) shows. In a neighbourhood of 
the surface in parameter space defined by Eq. (4.84) there exist supersymmetric minima 
with negative cosmological constant. We stress that the existence of these supersymmetric 
minima does not depend on the form of the Kahler potential and is, hence, not subject to 
the constraint (2.16) in moduli space. Consequently, these minima exist for an arbitrary 
number I of blow-up moduli and, in particular, for the specific model detailed in Table 1. 

In our discussion of the universal model with a single T and U modulus, we have encoun- 
tered three classes of supersymmetric minima characterized by the sign of the parameters. 
Above, we have shown that the analog of one of these classes, namely the one that includes 
supersymmetric minima with vanishing cosmological constant, also exists in the general case. 
What about the other two cases? We can construct examples for those cases by generalizing 
the results for the universal case. We choose specific parameters, such that 



ka — k 



ma 



m 



l = ll 
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7^' 



(4.86) 



for all a and i and universal parameters /c, 
the F-terms at the universal point 



m. 



I and /i. For these parameters, we compute 



T" = T , U' = U , 
in field space. We find that W = 7W where 

W = w{k,m,T) + w{l,fi,U) 

and 

Fa = w'{k,m,T)-^W 
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(4.87) 
(4.88) 

(4.89) 
(4.90) 
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where b = 81. Comparison with the universal F-terms (4.64) and (4.65) then shows that 
every solution (T, U) to the universal F-equations for parameters k, m, I, fi gives rise to 
a solution of the general F-equations via the identification (4.86), (4.87). We still have 
to show whether these solutions, if minima in the universal case, remain minima for the 
general model. To do this, we evaluate the expressions (2.22) and (2.23) for the second 
derivative of the potential at supersymmetric minima, using the general model defined by 
Eqs. (2.9) and (2.13) but specializing to universal parameters (4.86) and fields (4.87). The 
result is compared with that of an analogous calculation for the universal (T, U) model. This 
comparison shows that minima of the universal (T, U) model indeed remain minima of the 
general model at the universal point if t is sufficiently large and e small. As before, these 
minima still exist for non-universal parameters sufficiently close to the universal choice (4.86). 
We note that the phenomenological parameter b in the universal Kahler potential (4.59) is 
identified as 6 = 8/ under this correspondence, where we recall that / is the number of 
blow-up moduli. In our analysis of the universal case we found that consistent minima, 
stable under higher-order corrections to the Kahler potential, exist for b ^ 40 {b < 80 if one 
tolerates a small value of u). This translates into an upper bound of / ^ 5 (/ ^ 10 allowing 
small values of u) for the maximal number of blow-up moduli for which we can construct 
supersymmetric minima in this way. In summary, we conclude that the supersymmetric 
minima with negative cosmological constant can indeed be extended to complete models 
with up to 5 (or 10) blow-up moduli. 

Can the supersymmetry breaking minima with negative cosmological constant we found 
in the universal case be generalized to the full model? Let us consider a minimum for 
the universal model with parameters m, k, /i, / and field values T and U and analyze the 
general model at universal parameter values (4.86) and universal field values (4.87). It is 
straightforward to show that 

dV 

gp = 0. (4,91) 

using the fact that the U derivative of the universal potential (4.66) vanishes and identifying 
b = 81, as before. Unfortunately, the derivatives of V do not vanish exactly at the 
universal point, due to the non-trivial coupling between the T and U moduli encoded in the 
coefficients pia- However, it can be shown that dV/dT"- consists of terms either of order e or 
suppressed by inverse powers of t. Hence, for sufficiently large t and small e the derivatives 
dV/dT"- are small and there will be an extremum for field values close the universal choice. 
Moreover, under the same conditions on t and e these will still be minima. We conclude, 
that non-supersymmetric minima with negative cosmological constant can be generalized to 
the full model. As above, the constraint on the parameter b, necessary for consistent minima 
of the universal model to exist, translates into a bound of / ^ 5 (/ ^ 10 allowing small u 
values) on the number / of blow-up moduli for which we can obtain such non-supersymmetric 
minima. 
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5 Conclusion 



In this paper we have analyzed the vacuum structure of four- dimensional = 1 supergravity 
theories originating from M-theory on G2 spaces, with a superpotential from flux and mem- 
brane instanton effects. We have focused on G2 spaces which are constructed by blowing 
up G2 orbifolds. These spaces have two different types of moduli, namely "bulk" moduli T° 
associated with the underlying orbifold and "blow-up" moduli which measure the size 
and orientation of the blow-ups. 

We have been starting the analysis with a simple toy model consisting of a single bulk 
modulus T, generalizing to seven bulk moduli T", then including blow-up moduli, first in 
a universal model with a single bulk and blow-up modulus (T, f/), and, finally, studying 
the full model. Our main result is that minima with negative cosmological constant can be 
found. They exist for both supersymmetry preserved as well as broken and, after appropriate 
tuning of parameters, supersymmetric minima with vanishing cosmological constant exist for 
the general model, as well as for most of the simpler toy models. In constructing consistent 
minima we had to respect one technical constraints on the moduli space. The typical ratio 
e = u/t oi the a bulk modulus t and a blow-up modulus u had to be smaller than one since 
the Kahler potential has only been calculated to leading order in e. The supersymmetric 
minima with vanishing cosmological constant are unaffected by this constraint since they 
only depend on the superpotential. However, for all other minima it implies an upper bound 
on the number / of blow-up moduli. Moreover, we need a large parameters ka multiplying 
the Ta instanton contributions in the superpotential in order to generate sufficiently large 
values for the bulk moduli. We expect both restrictions could be avoided if the exact Kahler 
potential was known. 

Although we did obtain supersymmetry breaking minima we have not been able to find 
any examples with a positive cosmological constant. From our experience a positive cos- 
mological constant cannot be achieved in the given setting, using a combination of flux and 
membrane instantons. In analogy with the IIB construction of Ref. [21], this can presumably 
be achieved by adding wrapped (anti) M5-branes to our set-up. Work in this direction is in 
progress. 

Acknowledgments A. L. and B. d. C. are supported by PPARC Advanced Fellowships. 
S. M. is supported by a Royal Society postdoctoral fellowship. 

Appendix 

A Full potential for a single bulk and blow-up modulus 

We consider the case with a single bulk modulus T and a single blow-up modulus U, split 
into real and imaginary parts as 



T 



t + ir 
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U = u + iu . 

The model is defined by the Kahler potential and superpotential 

K = -71n(r + T)+6i^±^ 
W = mT + ke-^ + 12U + le-^ , 
where m, k, /i, / and b are real constants. We then find for the F-terms 
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The scalar potential is given by 
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